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We conduct extensive independent numerical experiments considering frictionless disks without 
internal degrees of freedom (rotation etc.) in two dimensions. We report here that for a large range 
of the packing fractions below random-close packing, all components of the stress tensor of wet 
granular materials remain finite in the limit of zero shear rate. This is direct evidence for a fluid-to- 
solid arrest transition. The offset value of the shear stress characterizes plastic deformation of the 
arrested state which corresponds to dynamic yield stress of the system. Based on an analytical line 
of argument, we propose that the mean number of capillary bridges per particle, v, follows a non¬ 
trivial dependence on the packing fraction, <f>, and the capillary energy, e. Most noticeably, we show 
that v is a generic and universal quantity which does not depend on the driving protocol. Using 
this universal quantity, we calculate the arrest stress, a a , analytically based on a balance of the 
energy injection rate due to the external force driving the flow and the dissipation rate accounting 
for the rupture of capillary bridges. The resulting prediction of o a is a non-linear function of the 
packing fraction f>, and the capillary energy e. This formula provides an excellent, parameter-free 
prediction of the numerical data. Corrections to the theory for small and large packing fractions are 
connected to the emergence of shear bands and of contributions to the stress from repulsive particle 
interactions, respectively. 

PACS numbers: 45. 70.Mg, 62.20.M-, 45.05.+X 
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I. INTRODUCTION 


Careful studies for dry granular materials suggest that 
the arrest of granular flow is universal in the sense that 
this non-equilibrium transition admits a continuum de¬ 
scription. This points to the existence of an equation 
of state for dry granular matter M- Similarly, vari¬ 
ous features of phase transition in vertically agitated wet 
granular materials can faithfully be described in terms of 
thermodynamic concepts [7H11|. 

Recently, a number of papers jT^ - ITTif also addressed 
phase transitions in sheared wet granular matter. In par¬ 
ticular, the arrest of wet granular flows driven by external 
forces with a cosine profile was attributed to a crossover 
of the power injected into the flow by the external field 
and the power dissipated in the rupturing of capillary 
bridges |l 2 l. IThi ITot . 

This provided a quantitative description of the arrest 
of shear flows of bidisperse disks [l 2 j and dumbbells [l 6 | 
with a fixed density, and of three-dimensional gravity- 
confined flow of monodisperse grains [l5|. The studies 
established that the minimal external forcing required to 
maintain flow can be calculated along the same line for 
transitions in fixed-density and fixed-pressure settings. 
Also, for both wet disks and dumbbells, a large hystere¬ 


sis was identified in which the fluid-to-solid arrest stress, 
cr Q , or the dynamic yield stress, has been found to be 
smaller than the solid-to-fluid yield stress, or the static 
yield stress. This is in accord with the conventional wis¬ 
dom where the stress at yield point is conceived to be 
larger than the minimum (plastic) stress which is re¬ 
quired to maintain the flow. Moreover, in a recent study, 
the difference of the static and dynamic yield stress is 
pinned down to the nonmonotonicity of the flow curves. 
Furthermore, it was shown that the nonmonotonicity of 
the flow curves indicates a shear banding instability 17 1 . 


Here, we augment these studies by addressing the ar¬ 
rest stress in systems with a prescribed global shear rate, 
7 . Rather than settings with complex flow profiles we 
consider for ourpresent study uniform shear flows in a 
Lees-Edwards [18|, [l9| periodic boundary flow geometry 
(Fig. [1]). This point of view is dual to the study of force- 
controlled systems where we prescribe a force amplitude 
and measure the minimal force required to sustain flow. 
In the present study we will address the arrest stress, 
i.e., the smallest value of the stress observed when de¬ 
creasing the shear shear rate in spatially uniform granu¬ 
lar flows. By definition of the stress tensor, the prod¬ 
uct of this arrest stress and the system size amounts 
to the minimal force that creeping flows with uniform 
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FIG. 1: (color online) Setup of the simulations with a snap¬ 
shot of the system. It comprises of a square box of length L 
sheared in vertical direction with a shear rate, 7 , and Lees- 
Edwards (skew periodic) boundary conditions. Disks are indi¬ 
cated by circles that represent their actual size. The strength 
of the forces, fij, acting between the disks i and j is indi¬ 
cated by lines connecting the disks, where the color marks 
the strength of the force according to the color bar to the 
left. On the color palette, negative (positive) numbers corre¬ 
spond to attractive (repulsive) force. The snapshot shows a 
uniformly sheared state of a simulation of size, L = 30, shear 
rate, 7 = 0.03, and capillary-bridge energy, e = 0.05, and 
packing density of disks, 4> = 0 . 8 . 


shear profile exert on the borders of a shear cell. This 
correspondence provides a quantitative prediction of the 
arrest stress where the only adjustable parameters have 
been measured in the force-controlled setting jl2|. The 
agreement provides further support for the modeling of 
the flow threshold based on energy dissipation arguments 
that were established in mill-' 

The present study is based on extensive simulations of 
a two dimensional wet granular material whose critical 
force for the arrest of flow has fully been characterized in 
Ref. [l2|. In contrast to that previous work we address 
now Lees-Edwards shear flow. We determine the depen¬ 
dence of the arrest stress, a a , on the packing density, <f>, 
and the capillary energy, e. 

The data are obtained by careful extrapolation of a se¬ 
quence of numerical measurements of the stress for larger 
shear rates to vanishing 7 , where we take care to base 
this extrapolation only on systems with uniform shear 
profiles. Requiring consistency with the setting where 
the flow is driven by an external force fl2l |. provides a 
parameter-free prediction of er a that is in excellent quan¬ 
titative agreement with the numerical results. 

The paper is organized as follows: In Sect. [TT] we re¬ 
visit the theory for the arrest of flow [l2j], and adapt it 
to describe the arrest stress of a uniformly sheared sys¬ 
tem. Subsequently, in Sect. cm we describe our system, 
its equations of motion, the approach adopted to solve 


them numerically, and the data analysis. Sect. IIVI com¬ 
prises the results on the shear-rate dependence of the 
stress, and its extrapolation to zero shear rate. We will 
demonstrate that for intermediate packing densities, <f>, 
the parameter dependence of the arrest stress is faith¬ 
fully described by the parameter free prediction derived 
in Sect. CD and trace down the additional physical pro¬ 
cesses causing the differences for small and large <fi. Our 
main results are summarized in Sect. m 


II. PREDICTING THE ARREST STRESS 


We consider uniform shear in a two-dimensional sys¬ 
tem in a domain of size L x L. By definition, the shear 
stress, a xy , corresponds to a force a xy L required to act 
at the boundaries of the system in order to maintain the 
uniform flow. Following [l2| we identify the critical force 
persisting at very small shear rates based on a power bal¬ 
ance of the work injected into, and dissipated in the sys¬ 
tem. Energy is dissipated by breaking capillary bridges 
spanning the stress network. Hence, the dissipated power 
takes the form 


(Pdiss} = / dxn s (x ) 

Jo 


dv y {x) 


dx 


ve = N 7 ve (1) 


where e is the energy needed to rupture a capillary bridge, 
and v the average number of bridges ruptures when two 
disks pass in the shear flow. To arrive at the flow con¬ 
figuration (i) the spatial density of disks, n s (x), per unit 
length is the ratio of the total number of disks, N, and 
the width, L, of the system, and (ii) the shear rate is 
constant, dv y (x)/dx = 7 =const. 

On the other hand, for a system with uniform density 
and shear the total injected power is the product of the 
off-diagonal element of the stress tensor, a xy , the shear 
rate, 7 , and the system area, Lr, 

(-^forcing) = &xy T A ■ (2) 


Flow ceases when the forcing injects too little energy to 
balance dissipation. The threshold value, a a , is obtained 
by balancing (Pdiss) and (Pforcing), 

cr a = n a ve (3) 

where n a = N/ L 2 is the areal density of disks. 

In this expression the number v depends on the pack¬ 
ing fraction, cj> , and on the capillary energy, e, because 
collective motion is needed when disks pass each other 
in a dense system—as initially discussed for a system 
where flow is driven bv an external cosinus-shaped force 
field [HI]. In Ref. [e|, based on systematic data col¬ 
lapses and heuristic scaling arguments, we have derived 
an explicit analytical formula for the dependence of the 
average bridge number, v, on the packing fraction, cf >, 
and the capillary energy, e. As an alternative equation, 
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FIG. 2: (color online) Fit of the amplitude, F, of the crit¬ 
ical cosinus force field required to maintain a flow in peri¬ 
odic boundary conditions. In Ref. ;12] it was shown that 
F = 8 ev/L, and we show here that the numerical data dis¬ 
cussed in that paper are well described by the expression, 
Eq. 0 , for v. The data shown by stars and filled circles refer 
to simulations at fixed e = 0.01 and packing densities <f> = 0.7 
and 0.82, respectively, where we varied the system size L as 
indicated by the top axis of the graph. The data shown by 
pluses and crosses show results of simulations for a systems 
size L = 18, capillary energies e = 0.01 and 0.05, respectively, 
and varying densities (f> as indicated on the bottom axis. Fi¬ 
nally, the open squares show data for L = 18, <j> = 0.7, and 
different values of e ranging between 10 -4 and 10 _1 . 
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in Fig. [2] we show that the data discussed in E2 are also 
well described by the following equation 



with B = 0.20 ± 0.02. 

In the following we explore in how far the prediction, 
Eq. ([31), complemented with the expression, Eq. Q, pro¬ 
vides a faithful description of the arrest stress in systems 
where the shear rate, 7 , is prescribed. 


III. NUMERICAL METHOD 

We numerically solve Newton’s equation of motion of 
a 1:1 mixture of large and small disks moving in a two- 
dimensional domain of size L x L subjected to Lees- 
Edwards boundary conditions El that enforce a pre¬ 
scribed total shear, L 7 , in the system. All velocities and 
forces are calculated based on the appropriate represen¬ 
tation of the flow in the two-dimensional plane that is ob¬ 
tained by properly unfolding the Lees-Edwards boundary 
conditions (c/. El). Crystalline ordering is prevented by 
considering a 1:1 mixture of large and small disks with a 
ratio Ri/R s = 1-4 of their respective radii Ri and R s [20j |. 
The random-close packing limit of this mixture has been 
reported |H |2l| to be <)) r cp = 0.84. We assume a con¬ 


stant mass per area, p, such that the mass of disk i is 
rrii = irpR 2 . 


A. Equations of motion 


We ado pt a fifth-order predictor-corrector Gear algo¬ 
rithm El El to solve Newton’s equation of motion 



] e ij ) 

ieAf(i) 


( 5 ) 


where rrii and 17 are the mass and the center disk i, re¬ 
spectively, A f{i) is the set of neighbors j interacting with 
i, the unit vector e, 7 points from the center of disk i to 
the center of disk j, and / ;J - is the force exerted by disk 
i on disk j. The latter force comprises a repulsion, and 
whenever applicable also the attractive force, F c . model¬ 
ing capillary bridges. 

Mutual repulsion between the disks is modeled by the 
repulsive force 


/1 / \ f Cij(Ri + Rj T ij ) for r ij Ri -f“ Rj 

r[nj) ~\0 else, 

( 6 a) 

where r ,; 7 is the Euclidean distance between he center of 
disk i and j. In the spirit of Hertz’s contact law El we 
set 


C i3 = C 


RiRj 

Ri + Rj 


1/2 


( 6 b) 


in order to account for different disk radii. The global 
parameter C controls the hardness of the disks. 

For the capillary bridges forces we engage the minimal 
capillary model proposed by Herminghaus in Ref. [24j 
in which the capillary force is constant, F c , as the dis¬ 
tance between granulates changes, and where it breaks 
at a distance s c . This model has been shown to be 
well-suited for numerical simulations of wet granular ma¬ 
terials ELI EE m. The force in this model is hys- 
teretic: upon first approach there is no attractive force 
(r'ij > Ri + Rj ); when the disks undergo collision (r,; 7 < 
Ri + Rj) a constant (capillary) force F c is activated, and 
the force persists until the disks separate by a (surface- 
to-surface) distance s c ( rij < Ri + Rj + s c ). At that 
point, T'ij = Ri + Rj + s c , the capillary bridge is re¬ 
moved (it ‘ruptures’). There is no force acting between 
the disks again, until they undergo their next collision 
{fij < Ri + Rj). Consequently, an energy e = F c s c is 
dissipated after rupture of the capillary bridge. 

Throughout this paper, we employ dimensionless 
rescaled quantities based on the capillary force, F c . the 
mass density of the disks, p, and the average disk di¬ 
ameter D. Time, t , and mass, m, is hence measured in 
units t = \J pD 3 /F Cl and p = pD 2 , respectively. Using 
these normalized quantities it is straightforward to nor¬ 
malize all physical quantities derived from mass, length 
and time, such as the local averages of disk velocities, 
components of the stress tensor, and the shear rate. 
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B. Data acquisition and analysis 


For each parameter set we run 20 simulations over a 
fixed total strain of 7 t = 10 x L. The choice to fix the 
total strain, rather than fixing the total simulation time, 
ensures that in the limit of very small flow rates we still 
sample a representative set of statistically uncorrelated 
disk configurations. We start to measure the physical 
quantities at the strain 7 = 4 x L and take 150 shots 
over the rest of the simulations until = 10 x L. Alto¬ 
gether these 150 configurations for 20 different runs pro¬ 
vide 3000 independent snapshots for each considered pa¬ 
rameter set. For each snapshot, we determine and store 
the components of the stress tensor, cr, as well as the av¬ 
erage density, partial densities, and granular temperature 
evaluated in ten bands parallel to the shear. 

The components of the stress tensor are calculated fl2l . 
mm by evaluating 


&a[3 — 


+ 


1 

1? 

1 

T? 


^ ' TTlj (Vi t a U a (pCi')') 
i 

^ ' r ij,oc Fij,f 3 • 
i<j 


K/3 - U P ( x 0) 


(7) 


In Eq. 0 a,/3 6 {x, y} denote the Cartesian components 
of the respective vectors or tensors. The stress tensor has 
two contributions: (I) the first term accounts for advec- 
tive momentum transport, where , is the velocity of 
the disk *, and U(x^) is the local drift velocity at the po¬ 
sition Xi of disk i, i.e., the overall center-of-mass velocity 
of the disks with centers in an interval of width L /10 
containing 2 j. (II) the second term describes the contri¬ 
bution of the interaction forces between particles on the 
stress to the flux of linear momentum. In this contribu¬ 
tion Tij is the vector connecting the centers of disk i and 
j , and Fij is the force exerted by the disk i on disk j. 
In the limit of small shear rates, 7 —> 0, one approaches 
creep flow, and the latter term in Eq. [7] always dominates 
the shear stresses. 

For non-Newtonian fluids, the relationship between the 
shear stress, a xy , and the shear rate, 7 , is given by 
Herschel-Bulkley relation [27} : 


(Txy = <7a + K X 7" (8) 

where, 07 , is the arrest (yield) stress, K, is a constant, 
and, n, is the shear thinning exponent. According to a 
recent review by Bonn et al. [28j . the shear thinning ex¬ 
ponent is a material dependent parameter for which no 
universality has been found. Therefore, Eq. [ 8 ] would be 
a natural choice to be considered as a template function 
for experimental data fitting. However, we found out 
that when the shear thinning exponent, n, is variable, 
for large capillary energies, e > 0.05, the fitting algo¬ 
rithm frequently takes very small exponents in the order 
of 10“ . This very small exponent gives rise to a flow 
profile which is very steep near the origin. As a result 
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FIG. 3: (color online) The upper panel shows the mean center- 
of-mass velocity in ten stripes parallel to the flow for three 
shear profiles encountered in simulations for shear rate 7 = 
0.02, system size L = 20, packing fraction cj> = 0.82, and 
capillary-bridge energy e = 0.1. The lower panel shows the 
respective relative deviations from the ideal profile. The open 
red squares connected by a solid line mark a profile where 
the deviations from the mean profile are in the admissible 
bounds of ±30% that are indicated by thin black lines. Green 
crosses and blue circles, that are connected by broken lines in 
according color, show profiles where the system has formed a 
shear band. 


of the sharp fall near the origin, it becomes more likely 
that the fitting algorithm grasps a negative offset. The 
negative offset corresponds to a negative arrest (yield) 
stress, cr a , which is non-physical. In order to prevent this 
fitting failure, it is required that the final fitting flow pro¬ 
file should be smooth at the origin. This is guaranteed by 
enforcing n > 1. To be consistent, we choose a quadratic 
profile, n = 2 , which is the most trivial choice given by 
the Bagnold scaling 53,121. 


C. Uniform shear profiles 

In the analysis of the arrest stress we systematically 
disregard all the snapshots where the velocity profile de¬ 
viates by more than a factor of 0.3 from the ideal linear 
case 

v y (x) =7 (x- , (9) 

hence excluding phase-separated states (c/. Fig. [3]). The 
ratio at which the snapshots are disregarded from further 
analysis strongly depends on the packing fraction, </>. The 
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higher the packing fraction, the lower the exclusion ra¬ 
tio. As an example, for (f> = 0.52, the rejection ratio 
is approximately ~ 72%, and for <f> = 0.83, it is about 
~3%. 

The lower panel of Fig. [3] shows the respective relative 
deviations from the ideal profile in order to exemplify 
the criterion for the selection of valid flow configurations. 
Flow configuration with deviations of the average flow 
velocities from the ideal profile (the thin dotted line in 
the upper panel of the figure) up to ±30% commonly 
appear as fluctuations on statistically uniform shear pro¬ 
files: these flow configurations will enter the analysis of 
the arrest stress. On the other hand, larger fluctuations 
rapidly evolve into shear bands that are persistent in 
time—once they appear, they do not decay again. The 
green crosses and blue circles, that are connected by bro¬ 
ken lines in according color in Fig. [3] show profiles where 
the system has formed a shear band. In that case the de¬ 
viations lie noticeably outside the admissible range, and 
they persist: the system has phase separated into two 
domains of almost constant drift velocity; localizing the 
shear in narrow bands that lie at 5 ~ x ~ 10 (green 
crosses) and 12 ~ x ~ 15 (blue circles), respectively. 
These configurations must be excluded from the analysis 
because they show a region with a very low (often zero) 
number of capillary bridges, and hence anomalously low 
shear stress. 


IV. RESULTS AND DISCUSSION 
A. Shear-rate dependence of the arrest stress 

In Fig. [fjlnset, we show the dependence of the com¬ 
ponents of the stress tensor on the shear rate j. The 
diagonal elements of the stress tensor, a xx and a yy , 
are depicted by red pluses and green crosses, respec¬ 
tively, and the blue stars show the off-diagonal element 
a xy = v y x- All curves show finite offsets at j = 0. 
The solid lines show the best fits by quadratic functions 
00 / 3 ( 7 ) = a a p + b a f 37 , where a al 3 accounts for the non¬ 
vanishing stress at zero shear rates, and b „jry 2 is the Bag- 
nold scaling expected to arise at large 7 [13, |29j] . In the 
main panel of Fig. [H we verify this quadratic scaling for 
large 7 by plotting the components of a — a on a double 
logarithmic scale. The offset of the off-diagonal element 
provides the arrest stress, <r 0 = a xy = a yx . It is an in¬ 
tensive parameter of the system that does not depend on 
the system size. However, clearly, the components of a 
and b are functions of the density, (j>, and the capillary 
energy, e. 

B. Parameter dependence of a a 

Fig. [5}a demonstrates our data collapse according to 
Eq. ([3]) for the dependence of the arrest stress, cr a , on the 
packing fraction, (f>, for five different capillary energies, e. 



shear rate 


FIG. 4: (Inset) The shear rate, 7 , dependence of the compo¬ 
nents of the stress tensor, a. The red pluses, green crosses, 
and blue stars show a xx , a yy , and a xy = a yx , respectively. 
The solid lines show the best fits by a quadratic function, 
o'api'y) = CLap+bapj 2 , to the data where, a xx — 0.095±277 2 , 
< T yy = 0.097 ± 277 2 , and a xy = 0.09 ± 67 s . One can see that 
all the components of the stress tensor possess finite offsets 
at 7 = 0. This refers to a finite dynamic yield stress of wet 
granular materials for (j> < cf) r c P ■ (Main Panel) Components of 
the stress tensor subtracted from their corresponding offsets, 
<7 — a a , in a double logarithmic scale. In these simulations, 
the packing fraction is <j> = 0.78, capillary energy is e = 0.05, 
and box length is L = 20. 


For 0.6 < (f> < 0.8 we find an excellent agreement between 
our data and the prediction cr a /(({>v£) = 1 shown by the 
dotted line. 

The data collapse is breaking down upon approaching 
the random close packing limit for values (j> > 0.8. We 
expect that in this range the repulsive forces gives rise to 
an additional contribution to the stress, as observed in 

0 • 

For cj) < 0.6 the values of 07 systematically fall be¬ 
low the prediction, Eq. ©■ This systematic error results 
from the choice a a ± b a j 2 of the function selected for ex¬ 
trapolating the 7 dependence of a xy to small 7 : it does 
not account for the slight decrease of a xy for small val¬ 
ues of 7 in the regime where one expects shear banding 
(c/. for instance Fig. 1 of 0 )- As shown in the main 
panel of Fig. []] there always is a small mismatch of the 
simulation data and the fit function for small 7 . We still 
adopted the function a + bj 2 because it involves fewer 
parameters and results in a more robust fitting of our 
data. To within our numerical error margins the system¬ 
atic error in the fit is negligible for </> > 0 . 6 , where the 
dip becomes so minute that shear banding is no longer 
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FIG. 5: (a) collapse of the data according to Eq. m, which 
asserts that (T a /{(t>vs) should be one (as indicated by the dot¬ 
ted line), as function of the packing fraction, (j>, for L = 20 
and five different capillary energies, e, that are provided in 
the lower left corner of the panel-b. (b) The panel provides 
the probability to observe shear bands in the simulations. In 
the legend the corresponding capillary energy, e, of each sim¬ 
ulation is given. 


encountered. On the other hand, for <p < 0.6 the fit 
provides a value characterizing the minimum of <r X y{i ) 
rather than the asymptotic value for vanishing 7 . When 
this distinction becomes noticeable the system becomes 
prone to shear banding ( Fig. ©b). 


V. CONCLUSIONS 

We calculated the arrest stress, a ai characterizing the 
minimal value of <r xy {j) for two-dimensional isotropic 
shear flows of wet granular disks by fitting the compo¬ 
nents of the stress tensor by a quadratic functions, a+ 67 2 . 
In the range 0.52 < <f> < 0.84, all components of the stress 
tensor of wet granular materials approach non-trivial fi¬ 
nite values for small 7 . The resulting offset value for the 
shear stress is of particular interest. In general it is differ¬ 
ent from the static yield stress, a y , that can be obtained 
M by considering the limit a y = lirmy_>o c xy . The latter 
characterizes the stress where a system at rest can start 
to flow. In the case of wet granular materials, the static 
yield stress is the force which is needed to break down 
the percolation network of capillary bridges. The former, 
i.e., the dynamic yield stress, is the minimal value of the 
stress which is required to maintain the flow. 

The prediction of the arrest stress is based on a bal¬ 


ance of the power injected by increasing strain at a fi¬ 
nite stress and the energy dissipation due to rupturing 
capillary bridges. For uniform shear this argument pro¬ 
vides Eq. ©. This prediction has a single free parameter, 
the average number, v, of bridges ruptures when parti¬ 
cles pass each other in the flow. We have determined 
this number fl 2 | in the context of a different flow, non- 
uniform shear flow driven by an external force field in 
periodic boundary conditions. Evoking universality of 
the number, v, with respect to the cause of the flow pro¬ 
vides a prediction of the arrest stress without adjustable 
parameters. The excellent data collapse shown in Fig. © 
a confirms the versatility of this approach to determine 
the arrest of flow [nmm, and the universality of v as 
far the protocol of inducing the flow are concerned. The 
difference between our prediction, Eq. m, and the nu¬ 
merical findings for <j> > 0.8 have been associated to the 
emergence of contributions to the stress from repulsive 
forces when approaching random close packing. The sys¬ 
tematic difference for </> < 0.6 is due to a noticeable devi¬ 
ation of the arrest stress from the static yield stress. The 
larger the difference, the more likely the flow will phase 
separate into a state featuring shear bands (c/. Fig.©b). 

Previous studies of far-from equilibrium transitions in 
wet granular materials have unraveled origins of these 
transitions to be either force- or energy-driven I EH- For 
instance, for wet granulates in a Petri dish under vertical 
vibration, solid-to-fluid (fluid-to-gas) transition is force- 
driven (energy-driven). In sheared wet granular matter, 
the solid-to-fluid transition, whose threshold represents 
the static yield stress, has been shown to be of force- 
driven origin [30] 

The present work complements and further substanti¬ 
ates the findings of previous studies nnnm that at¬ 
tributed the fluid-to-solid transition in sheared wet gran¬ 
ular systems is connected to a balance of the power in¬ 
jected by the forcing of the flow and the energy dissipa¬ 
tion rate due to rupturing of capillary bridges in a plastic 
shear flow (energy-driven). This mechanism is funda¬ 
mentally connected to the fact that wet granular mate¬ 
rials have an inherent energy scale, the energy e that is 
dissipated upon rupturing a capillary bridge 7.E3- it 
will be interesting to explore the impact of this impor¬ 
tant distinction for other soft matter systems [ 13 , |3ll 0 ] 
where small attractive interactions and a dissipation due 
to transfer of kinetic energy into internal degrees of free¬ 
dom of the system should admit a corresponding analysis 
of the arrest of flow. 

This is in contrast to dry granular materials where the 
arrest of flow is expected to be controlled by frictional 
forces between the particles. 
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